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We present the Schmidt decomposition for arbitrary wave- 
functions of two indistinguishable bosons, extending the re- 
cent studies of entanglement or quantum correlations for two 
fermion systems [ J. Schliemann et a i, Phys. Rev. B 63, 



085311 (2001) and |quant-ph/0012094| . We point out that the 
von Neumann entropy of the reduced single particle density 
matrix remains to be a good entanglement measure for two 
identical particles. 
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Recently, quantum information and computing has 
emerged as an active research field. Of the more fun- 
damental issues in quantum information science, charac- 
terization of inseparable quantum correlations, or entan- 
glement has received considerable attention . As is well 
accepted by now, the power of quantum computing lies 
mostly at the inseparable nature of multi-partite quan- 
tum states. In fact, it is often said that entanglement 
is a quantum information processing resource (currency) 
which all known quantum information protocols need to 
consume. Historically, quantum entanglement was first 
raised by Schrodinger |Q in responses to EPR's local re- 
alism Despite more than 70 years of successful ap- 
plications of quantum theory, however, we still do not 
understand the full content of this quantum inseparabil- 
ity. At the simplest level of multi-partite pure states, a 
complete knowledge for classifying and measuring entan- 
glement exists only if two distinguishable particles are 
involved. For mixtures or states of more than two par- 
ticles, our current understanding of inseparable correla- 
tions remains limited 

For quantum computing applications, generation of en- 
tanglement is the prelude of any useful protocol [Q . One 
commonly adopted mechanism relies on a direct interac- 
tion of the constituent quantum particles For indis- 
tinguishable particles, e.g. qubits represented by bosonic 
or fermionic atoms, photons, or electrons as often envi- 
sioned in scalable quantum computer architecture a 
direct short-range interaction may cause significant over- 
lap of individual particle wavefunctions, thus compro- 
mising the ability of keeping track of which particle is 
which qubit. This calls for understanding inseparable 
quantum correlations for identical particles including ex- 
change symmetry. 

The aim of this study is to investigate entanglement 
measures for systems of two identical bosons. We ex- 
tend the recent results on two fermion systems by J. 
Schliemann et al. Wpl, and present a unified entangle- 



ment measure in terms of the von Neumann entropy of 
the reduced single particle density matrix. This paper 
is organized as follows: we first present the well-known 
result for two distinguishable particles. We then summa- 
rize the recent results of two fermions and show that 
the von Neumann entropy of the reduced single particle 
state remains a good entanglement measure. Finally we 
provide a simple proof of the Schmidt decomposition of 
two boson wavefunctions and obtain its von Neumann en- 
tropy measure of entanglement in terms of the reduced 
single particle density matrix. We conclude that in all 
three cases: two distinguishable particles, two fermions 
or bosons, the reduced single particle density matrix fully 
characterizes their inseparable quantum correlations. 

Physically, we arrive at the following conclusion: for 
two identical particles, being simply in usually insepara- 
ble quantum states does no necessarily imply their use for 
quantum information science. In agreement with the au- 
thors of Ref . 1^ , we will use quantum correlations, rather 
than entanglement, to characterize useful inseparable cor- 
relations of two identical particles that are beyond what 
is normally required by the permutation symmetry. 

Two distinguishable particles For two distinguish- 
able particles A and B, any joint pure state wavefunction 
can be expressed as 
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where C is a complex matrix (NxM). aj (a^) and 5] {bj) 
are respectively the single particle creation (annihilation) 
operators for A and B in orthogonormal single states 
\i)A (of dimension N) and \j)B (of dimension M). In the 
position representation, one simply has aj|0) — + ipii'^A) 
and ^]|0) — *■ 4>j{rB)- The normalization condition gives 
tr(CtC) = tr(CCt) = 1. 

The Schmidt decomposition of (|l|) can now be con- 
structed pHl0|]. Using the theorem of Singular Value 
Decomposition (SVD) for a complex matrix we 
write C — U^DV with unitary matrices U {N x N) 
and V {MxM), and a diagonal matrix D [of dimen- 
sion min(A^, M)]. Transforming the single particle basis 
according to 

N 

1=1 

M 
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the state vector (|l|) becomes diagonal 

min(A',A/) 

Yl D,a2b2\0). (2) 
fc=i 

Therefore, it is said to be entangled if more than one 
diagonal elements of D (Dk) is nonzero, i.e. if the rank of 
matrix C is larger than 1. A good entanglement measure 
in this case is the von Neumann entropy of the reduced 
one particle density matrix. The two particle density 
matrix is p — and the reduced density matrix 

is defined according to = tTsip) and p^ = tr a{p), 
whose elements can be evaluated according to 

= i^HKl^) - (C^C),,.. (3) 

We immediately note that both Schmidt coefficients 
and the rank of matrix C are independent of uni- 
tary single particle transformations on I^P). Furthermore, 
the von Neumann entropy measure for entanglement 

S = -trA[/3Aln(/3A)] = -trs[/Os InCps)] 

^- Yl (4) 

is also representation independent. Because < |D^| < 1 
and J2k \^k\ — 1' is ^^^^ non-entangled states and 
takes a maximum value ln[min(7V, M)] for the maximum 
entangled state. 

Two-fermions A general two fcrmion state can be 
expressed as 

N 

1*^) = Y ^^,fIfh^). (5) 

where w is a complex and anti-symmetric matrix (NxN) 
LJij = —ujji. fl (fi) is the creation (annihilation) 
operator of a fermion in single particle state \i), and 
in position representation ///J|0) [V'i('^4)V'j(^s) — 
V'i(rB)'0i(rA)]/'\/2. The normalization condition is 
tr(wtt^) = 1/2. 

As was first shown in Ref . |7| , and extended in Ref. [^j , 
a Schmidt decomposition also exists in this case (^). A 
unitary transformation of fi — ^ijf'j ^^"^ fermionic 
operators f'^ and /j, transforms w in a non similar way. 
The state (||) becomes 

\^F) = f2 (c^w*)^^./;t/;io). (6) 

There exists a matrix U such that uj takes a block diag- 
onal form , containing 2x2 blocks of the type 



with eigenvalues Zk- Consequently the state becomes 
'diagonal' § 

< — 

\^F)=2Yzkf!,i_j;,l\0). (8) 

fc=i 

Several earlier studies on two electron correlations in 
atomic systems also used this Slater-Schmidt expansion 
]l3| , p^ . In Ref. Grobe et al. called it the 'canonical 
representation' (see also jl^), and a correlation measure 
was defined to he K = l/ix^pf^). Approximately, K 
measures the number of Slater orbitals needed for the 
expansion of a two electron wavefunction. K ^ \, there- 
fore, corresponds to no quantum correlation. 

Schliemann et al. ||^,^ propose to use the Slater rank, 
the number of non-zero z^'s, as the criterion for quan- 
tum correlations. If the Slater rank is larger than 1, 
the two fermion state 1 5* p) contains additional quantum 
correlation beyond the anti-symmetric permutation re- 
quirement. Immediately, one sees that there is no quan- 
tum correlations for two fermions in two or three states 
iV = 2, 3. The simplest case where additional inseparable 
quantum correlations can exist occurs for A'^ = 4. Sev- 
eral entanglement measures were proposed and studied 
in Refs. for this case. 

Compared to states for two distinguishable particles, 
we note the analog of the von Neumann entropy remains 
a good correlation measure for two fermions p5[ |. The 
two fermion density matrix is pp = |^'f)(^_f|- Because 
the two fermions are identical particles, there is now only 
one single particle density matrix, p^ (normalized to 1), 
which can be computed according to 

^l{^F\flU\^F)^2iJiu)^,. (9) 

The normalization condition is thus reduced to simply 

'l2k=i'^ kitP = 1/4- The corresponding von Neumann 
entropy is then 

^/=-tr[p^ln(p/)] 
<f 

= -ln2-4^|zfe|2ln(|zfen. (10) 
fe=i 

It provides a smooth measure from the non-correlated 
state Sf = In 2 to the maximum correlated state with 
Sf — Iu^Ne)- Ne is the largest even number not larger 
than N, i.e. Ne ^ N for even N and A^^; = - 1 for 
odd A^. 

We note that both the Slater rank and the von Neu- 
mann entropy Sf are independent of unitary transforma- 
tions on the single particle basis. 

Two bosons We now study wavefunctions for two in- 
distinguishable bosons 
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In this case (3 a complex symmetric matrix (NxN) j3ij = 
Pji- bl (bi) are the boson creation (annihilation) oper- 
ators for single particle state \i), and in position repre- 
sentation blb'^j_^.\0) [ipt{rAHjirB) + '4'i{rB)tpj{rA)]/V^ 

and 6|5||0) — *■ '!/'i('Vi)V'i(^s)- The normalization gives 
tr(/3^/3) = 1/2. Under a unitary transformation of &'s to 
new bosonic operators 6"s according to bi = Uijb'p 
the matrix (3 undergoes a non-similarity transformation 
(3 = W(3U* . The new state vector becomes 



N 



We now prove that the Schmidt decomposition also exists 
in this case. According to SVD , a complex symmetric 
matrix needs only one unitary matrix U , such that /3 = 
UBU^, and B being a diagonal matrix (with complex 
diagonal elements Bk). The proof is rather simple. SDV 
gives P = UBV"^, while /3 being symmetric gives f3 = 
0^ = VBU'^. Since SDV is unique, we obtain U = V 
I p^ . In fact, the matrix U may be chosen to be the one 
that diagonalizes j3. Thus the state vector of two bosons 
can be represented in a diagonal form 



N 



\^B)=V2Y,Bub'lb2\Q). 



(13) 
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We note that such a Schmidt expansion was suggested 
earlier in Refs. ||l^,|lj]. The inseparable correlations cri- 
terion is again related to the number of non-zero diagonal 
elements i.e. the rank of matrix (3. In general quan- 
tum correlations can arise as long as iV > 2. 

Again there is only one single particle density matrix 



{^B\blbA^B) 



= 2(/3^/3)^.. 



(14) 



Similarly to the fermionic case, the von Neumann entropy 



N 



Sb - ^tripHnp"] = -^(2|Bfc|2)ln(2|Bfe| 



(15) 



fc=i 



measures quantum correlations between the non- 
correlated state Sb — and the maximum correlated 
state Sb = In(iV). 

Finally we note as discussed in Refs. within a 

given rank, the determinant of the corresponding non- 
zero sub-matrix can also be used to measure quantum 
correlations. This was first shown in the simplest case of 
iV = 4 for two fermions (7| by introducing a dual space 
(dual matrix for lo). One can show in general the sub- 
matrix determinant takes the maximum value when Sf 



or Sb is the maximum. We provide the simple case of 
two bosons below. Let's assume the rank of /? matrix rb 
is equal to N. We denote \Bh\'^ — and construct a 

fimctional = 0/0=1 — HJ2k=i — 1/2)- The second 
term is due to the requirement for a proper normalization 
of the matrix /3^/3. We find the condition for a maximum 
by taking partial derivatives d^-J- — 0. The resulting 
equations then become Y[i=ij d — ^- Therefore the ex- 
treme occurs at £.1 = £.2 ~ ■ ■ ■ = £.N = 1/2 A^. We may 
then express the determinant of f3 as 



N 



^=ei(6,6,...,^jv)n^^- 

J=2 



(16) 



(12) Its second order derivative at £i — 1/2N is then 



- (% - 2) 



1 

27V 



N-l 



(17) 



which forms strictly negative matrix. Therefore the sub- 
matrix determinant does take its maximum at this ex- 
treme point of S*;,, a point which corresponds to the max- 
imally correlated state within the same rank. Extending 
the proof to rb < N is straightforward. 

Finally, we note that for two bosons our proposed en- 
tanglement measure in terms of the single particle den- 
sity matrix shares same spirit with the so-called con- 
densate fragmentation phenomena in systems of many 
bosons . Therefore it may also provide a useful start- 
ing point for the investigation of quantum correlations 
beyond two particle systems. One obvious consequence 
for many-boson systems is that fragmented condensates 
contain non-zero inseparable quantum correlations or en- 
tanglement. 

From a technical point of view, we note Schmidt de- 
compositions for two particle pure states can be con- 
structed starting from unitary transformations that di- 
agonalize the respective reduced single particle density 
matrices. We hope this unified approach will stimulate 
further discussions into quantum information science of 
indistinguishable particles. In the next paper, we will 
address the question of such inseparable quantum corre- 
lations for two bosons in quantum mixture states. 

In conclusion, we have presented the Schmidt decom- 
position for pure states of two indistinguishable bosons. 
We find the rank of the Fock space wavefunction matrix 
can be used to determine if additional quantum corre- 
lations (or entanglement) beyond the usual permutation 
symmetry exists. We have suggested the von Neumann 
entropy for the reduced single particle density matrix as a 
smooth entanglement measure for both two fermion and 
two boson wavefunctions, just as in the case of two distin- 
guishable particles. There remains one point of concern 
in this von Neumann entropy measure in that it did not 
seem to reflect the information content of quantum cor- 
relations. For instance, the un-correlated two fermion 
state gives Sf = In 2, rather than as in the classical 
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case. This points to the need for studies of quantum in- 
formation coding and processing with correlated identical 
particle states. 
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Note added. — After this paper was submitted for re- 
view, a preprint on the same topic of entanglement 
in a two-identical-particle system appeared. It provided 
a different proof of Eq. (|lj) in agalogy with the orig- 
inal approach used for two fcrmions |^]. The reduced 
single particle density matrix was not invoked in p7| . 
Although both results are essentially the same in terms 
of mathematics, we seem to arrive at different conclu- 
sions. In particular, an interesting construction is pro- 
vided in Ref. jlTt which transforms a symmetric state 
\Bk\{e"l"'b^^b^^ + e-"t"'h^^h^^)\0) of a degenerate pah 
of non-zero diagonal element Bk ~ \Bk\e'^'^^ into state 
^Ii^fc2l0)- Since both states are of rank 2, we consider 
them as containing quantum correlations. This is further 
evidenced by the fact that both can be used for quantum 
teleportation ||l8|. On the other hand, Ref. considers 
them as separable according to their definition 1. 
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